Nonlinearity exponents in lightly doped Conducting Polymers 
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The I-V characteristics of four conducting polymer systems like doped polypyrrole (PPy), poly 
3,4 ethylene dioxythiophene (PEDOT), polydiacetylene (PDA) and polyaniline (PA) in as many 
physical forms have been investigated at different temperatures, quenched disorder and magnetic 
fields. Transport data clearly show the existence of a single electric field scale for all systems. Based 
upon this observation, a phenomenological scaling analysis is performed, leading to extraction of 
a numerical value for a nonlinearity exponent called xm which serves to characterize a set of I-V 
curves. The conductivity starts deviating from an Ohmic value ctq above an onset electric field F 
which scales according to F a ~ <7q m . The electric field-dependent data are shown to be described by 
Glatzman-Matveev multi-step tunneling model [JETP 67, 1276 (1988)] in a near-perfect manner over 
nine orders of magnitude in conductivity and five order of magnitudes in electric field. Furthermore, 
xm is found to possess both positive and negative values lying between -1/2 and 3/4. There is no 
theory at present for this exponent. Some issues concerning applicability of the Glatzman-Matveev 
model are discussed. 



PACS numbers: 72.20.Ht,72.80.Le,05.70. Jk 



I. INTRODUCTION 



Conducting polymers (CPs) such as doped poly- 
acetylene (PA), polypyrrole (PPy), poly 3,4 ethylene 
dioxythiophene (PEDOT) and polyaniline(PANI) show a 
great variety of transport properties. In general, the elec- 
tronic structure of 7r-conjugated pristine (undoped) con- 
ducting polymers originates from the sp 2 p z -hybridizcd 
wave functions of the carbon atoms in the repeat unit. 
Despite strong electron-phonon coupling in pristine CPs, 
an extraordinary large range of conductivities has been 
covered by doping. Conductivity a range from the highly 
insulating values 10 -6 S/cm to highly metallic ones 10 5 
S/cm depending on the doping concentration 1 . Such a 
wide range of conductivities have made CPs useful in 
many applications such as wires, electromagnetic inter- 
ference shields, anti-static coatings, conducting layers in 
active devices including organic- and polymer-based light 
emitting devices, photovoltaic devices, and field-effect 
transistors^—. 

The electronic transport properties of CPs are fur- 
thermore strongly influenced by the synthesis procedure 
and intrinsic disorder^. The conductivity of CPs origi- 
nates from mobile charge carriers into the 7r-conjugated 
electronic orbitals which get filled by doping. At low 
doping densities, these charges self-localize to form soli- 
tons, polarons and bipolarons 7 . At higher doping levels, 
a transition to a metallic state is typically observed 8 -. 
Quite often, disorder dominates the macroscopic proper- 
ties, thereby hiding or even eliminating the intrinsic de- 
localization along the chains^. As developed in Ref. ITol . 
the intrachain transfer integrals are large as compared 
to the interchain transfer integrals, such that the trans- 
port takes place preferably along the chains. However, 
no matter how well the chains are ordered, the electrons 



finally have to hop between chains, and the related lo- 
cal resistances strongly influence or even dominate the 
macroscopic conduction. This leads to an increase in 
conductivity as the temperature is increased, i.e., an in- 
sulating behavior. Metallic behavior is observed only at 
larger temperatures and only in some polymer systems^. 
Thermoelectric power(TEP) measurements as a function 
of temperature, as well as anomalous differences between 
optical and d.c. conductivity data, provide evidence 12 
that CPs can be treated as being structurally hetero- 
geneous, consisting of thin, metal-like fibrils, separated 
by potential barriers^. This structural disorder is re- 
sponsible for low electrical conductivities owing to low 
mobilities for most conducting polymers, even though 
in fully doped polymers, charge carrier concentrations 
can be as large as 10 23 /cm 3 , which is about four or- 
ders of magnitude higher than in highly doped inorganic 
semiconductors. However, highly doped CPs show an 
intrinsic metallic nature with traditional signatures like 
a Drude metallic respons o 13 ' 14 , temperature independent 
Pauli susceptibility 1 — and a linear dependence of the ther- 
moelectric power on temperature 16 . For the majority of 
CPs, the metallic state is strongly dependent on struc- 
tural disorder, and only a small fraction of carriers re- 
main delocalized at low temperatures with long scatter- 
ing times. These delocalized electrons generate a weak 
temperature dependence of a down to millikelvin tem- 
peratures. 

Polypyrrole has been particularly extensively investi- 
gated due to its many prospects for applications, like 
relatively high environmental stability, high conductiv- 
ity, or the simplicity of preparation either by chemical 
or by electrochemical polymerizatio n 17 ' 18 . Polypyrrole is 
an amorphous conjugated polymer based on an aromatic 
ring and has a nondegenerate ground stated. The poly- 
mer chains are intertwined and the fibrils are randomly 
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oriented. Consequently, polypyrrole systems must be re- 
garded as three-dimensional disordered systems with re- 
spect to their structure and morphology. Over the last 
few years, PEDOT has attracted a lot of interest because 
of high conductivity, optical transparency, easy process- 
ability and high stability 20 . Inganas and co-workers 2 ^ 
have shown that PEDOT has a band gap of approxi- 
mately 1.6 eV and can be cycled between the reduced 
and the oxidized state. It has been suggested that PE- 
DOT has a lamellar-type structure built from ellipsoidal, 
conductive particles 22 . Even in highly conducting sam- 
ples, one observes a mixture of insulating and metallic 
transport behavior. This has been explained in terms of 
a heterogeneous morphology 1 ^. 

Conducting polymers are known to undergo a metal- 
insulator transition as a function of doping 23 . Samples at 
the insulating side of the transition are of interest here. 
Conduction data^ 4 - - — in CPs at low fields in this regime 
are usually discussed within the framework of the stan- 
dard Mott variable-range hopping (M-VRH) model— and 
its modification (ES-VRH) by Efros and Shklovskii 3 - - in 
presence of Coulomb interactions at low temperatures. In 
any VRH model, an electron near the Fermi level hops 
between localized states by absorbing phonons and the 
conductivity er is assumed to follow from path optimiza- 
tion of hopping distance and difference in energy levels 
involved: 



cr(T) = cr m exp 



(1) 



Here a m is the conductivity pre-factor and T the tem- 
perature. In the M-VRH model, the exponent m is given 
by m = l/(d + 1) where d is the dimensionality. Thus, 
m is 1/4, 1/3 and 1/2 in three, two and one dimension 
respectively. The characteristic temperature T a is given 
by T = 4Tr/3N(E F )k B a 3 where N(E F ) is the density of 
states at the Fermi level, ks the Boltzman constant and 
a the localization length. In the ES-VRH model, interac- 
tions among electrons were shown to open up a coulomb 
gap in the density of states at the Fermi level. This leads 
to m = 1 /2 independent of the dimension and 



T = 



2.8e 2 



(2) 



where e is the dielectric constant. The standard VRH 
model seems to describe well three dimensional systems 
like powdered FeCl 3 -doped PPy 2 ^, PEDOT at higher 
temperature 2 ^ yielding m=l/4. PTS-doped PPy films 26 
are reported to yield m=l/3 suggesting two-dimensional 
nature. The exponent in single crystals of PTS-doped 
polydiacetylene (PDA)2£ is found to be ~0. 65-0. 70 ap- 
parently demonstrating quasi-one dimensional nature of 
the system. On the other hand, some samples like PE- 
DOT at low temperature 2 ^, PPy films 2 ^ tend to follow 
ES-VRH particularly at low temperature with m=l/2. 
Apart from the VRH model, another model 31 called 
fluctuation-induced tunneling (FIT) model also has been 



used to describe electrical conduction in CP's, especially 
in quasi one-dimensional structures such as nanotubes or 
nanofibres. This model characterizes electrons transfer 
across the insulating barriers between conducting regions 
and hence, can be appropriate for describing the trans- 
port properties of CPs as individual polymer fibrils have 
a mixture of both crystalline and non-crystalline insulat- 
ing regions. In this model, the temperature dependence 
of a is given by 



cr(T) = crmexp 



T s + T 



(3) 



where is a temperature-independent prefactor, T rep- 
resents the typical magnitude of tunneling barriers, and 
the ratio T s /T determines the extent of reduction in con- 
ductivity at low temperatures. The model, originally de- 
veloped to explain the temperature dependence of con- 
ductance in C-PVC composites^ 1 -, was subsequently ap- 
plied to doped polyacetylene 3 - 2 -, isolated nanotubes or 
nanofibre o 33 ' 34 . It is known that the latter often exhibit 
conduction properties different from those of its assem- 
blies. The two models differentiate themselves in a fun- 
damental manner by their behavior as T approaches zero 
- the VRH conductivity goes to zero even as the FIT one 
remains finite. As seen later, this will have a bearing on 
the nonlinear behavior. Note that both the models in 
special cases (i.e. m=l in VRH or T s = in FIT) can 
represent an activated process. In rest of the paper, we 
will use the symbol a for conductivity and £ for con- 
ductance. The latter here is simply taken as the chordal 
conductance, I/V unless otherwise mentioned. 

Non-Ohmic or nonlinear conduction is a very common 
feature of disordered systems in general, and conducting 
polymers in particular. Onset of nonlinearity in these 
systems often takes place upon application of only few 
volts of bias across samples in laboratories. In bulk sys- 
tems, this is rendered possible by the presence of mi- 
croscopic inhomogeneties which lead to very large local 
fields. Now, with advent of low dimensional materials 
such as nanotubes, nanofibres, nanowires or quantum 
dots, the generation of large fields with moderate bias 
applied across small lengths are quite common. This pro- 
vides added incentive to study and understand the phe- 
nomena of nonlinear conduction in CPs. Nonlinear trans- 
port data in conducting polymer o 18 ' 27 ' 28 ' 34 " — are either 
presented in form of I-V curves or equivalently, in form 
of E-V curves, measured as a function of temperature. 
However, there is hardly any systematic analyses of the 
data mainly because theoretical understandin g 31 ^ 5 ' 40 " — 
is far from complete. For example, the models of vari- 
able range hopping under fielder— have expressions for 
the field-dependent conductivity in two field limits - low 
or moderate and high - as a function of temperature. 
Glazman and Matveev (GM)^ 4 . have evaluated a model 
involving multi-step tunneling across thin disordered re- 
gions also in two limits, eV 3> ksT and eV -C ksT . In 
the first limit, one has at any temperature a full expres- 
sion for E(V) containing parameters without having any 
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temperature dependence. 

In this paper, we present electrical transport mea- 
surements covering both linear and non-linear regimes 
in moderately doped PPy in forms of pellets as well 
as films and in a pellet of doped PEDOT over a wide 
range of temperatures. In addition, some published data 
from literature have been also processed for compari- 
son. For analysis of our data, we depart from methods 
used in existing literatures and adopt a newly suggested 
methodology^ based upon observation of existence of a 
field scale in many disordered samples. This method of 
scaling analysis of a set of I-V curves yields an expo- 
nent xt, called nonlinearity exponent, that characterizes 
the nonlinear data and is believed to reflect the under- 
lying field-dependent conduction mechanism. We apply 
this scaling approach to various conducting polymer sys- 
tems, examine its validity and extract the nonlinearity 
exponents. In the next section (section II), the method 
of scaling analysis is described in detail. Furthermore, 
two models which are extensions of those earlier consid- 
ered at low fields, namely, VRH and FIT models are also 
reviewed along with corresponding experimental results 
reported in the literature. The GM-model^ 4 - was origi- 
nally developed for thin amorphous films whose dimen- 
sions along the hopping direction (i.e. thickness) lie be- 
tween the localization lengths and hopping lengths and 
has never been used before to explain nonlinear data in 
CPs. It is found rather surprisingly to be very success- 
ful in explaining the I-V characteristics in CPs. These 
models are critically discussed in particular context of 
the general requirements of scaling. Experimental details 
are given in section III. Data are presented and analyzed 
mainly by adapting the GM-expression as a scaling func- 
tion in section IV. Results so obtained are then discussed 
in section V. Finally, conclusions are given in section VI. 



II. SURVEY OF MODELS OF NONLINEAR 
CONDUCTION: THEORETICAL AND 
PHENOMENOLOGICAL 

As the form of disorder can be of huge variety, so could 
be the mechanisms causing deviations from Ohmic be- 
havior. Here we will concern ourselves with the type of 
non-Ohmic phenomena that arise in bulk due to inter- 
play of disorder and various fundamental processes such 
as localization, tunneling, hopping etc. Thus, we exclude 
from our purview cases like space charge related conduc- 
tion in organic solids and joule heating, that also occur 
in bulk. The joule heating is, of course, present in any 
conductor, whether ordered or disordered, and increases 
with passing current. In case of positive temperature co- 
efficient of resistance (as in metals), the Joule heating 
leads to an increase in resistance whereas in case of neg- 
ative temperature coefficient of resistance (as in systems 
under consideration), it leads to a decrease in resistance. 
Since both the applied bias and Joule heating can lead to 
decrease in sample resistance, care must be taken during 



experiments to ensure that effects of joule heating is in- 
significant or taken it into account in the interpretation 
of data. 

A study of field-dependent conduction is generally ex- 
pected to bring out subtleties or additional processes 
that are either not present or insignificant in the lin- 
ear conduction. Consider, for example, composites which 
are random mixtures of conductors and insulators above 
the percolation threshold. For small bias, conduction 
takes place only through the backbone and is Ohmic. As 
the bias is increased, tunneling across the thin insulat- 
ing layer between the ends of the conductor chains dan- 
gling off the backbone are believed to lead to increased 
conductance 4 ^. Furthermore, a theoretical treatment of 
non-ohmic conduction may lead to self-consistent rela- 
tions among parameters that are also involving in Ohmic 
conduction. This is illustrated in case of VRH under field 
(see below) where physical parameters like localization 
length may be involved in more than one experimentally 
measurable quantities. In view of the inadequate theo- 
retical development it should be of much utility to know 
general properties of nonlinear conduction. One can then 
check whether any specific theoretical model prediction is 
consistent with general requirements or not. One general 
property is that the conductance £ always increases with 
field F, at least at small fields. This property is obvious 
since the application of a field results in lowering barrier 
heights, thus decreasing resistance to conduction. For a 
proper description, one defines a characteristic field or a 
field scale F such that > S(0) for F > F a . This 

corresponds to the fact that a sample remains Ohmic at 
small fields and starts deviating from the Ohmic behavior 
as the field is increased beyond F a . The scale may also be 
given formally by the following (Acr(F) = u{F) — c(0)): 

Aa(F ) ~ ,7(0). (4) 

I- V curves are generally measured as a function of some 
physical parameter upon which er(0) depends. Tempera- 
ture is the most commonly used parameter. Others like 
magnetic field, pressure are also valid parameters. An 
interesting but less studied situation is when I-V data 
are gathered simply as a function of quenched disorder 
as in composite samples 45 . Naturally, F Q is expected to 
be a function of the parameter. 

1. Scaling analysis: Nonlinearity exponent 

Let us consider, for definiteness, cases where measure- 
ments are carried out at different temperatures. Re- 
cently, a generalized approach based upon existence of 
a single field scale in a given disordered sample has been 
adopted, leading to characterizing a set of I-V curves by 
a single number xt, called nonlinearity exponent^. In 
this approach, the conductivity cr(T, F) is given by the 
scaling relation: 



<r(T,0) 
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where g is a scaling function and the field scale F a (T) 
at each temperature is given by the phenomenological 
relation 

F {T) = A T a Q XT , (6) 

where u (T) = cr(T, 0) is the conductivity at zero bias 
at temperature T and At is a constant whose value de- 
pends upon the criterion that fixes the scale F a . For 
F < F Q , the scaling function g s» 1 corresponding to the 
fact that the conductance increases very little from the 
zero field value a . At larger fields F > F a , g > 1. Thus, 
the field F Q can be called a crossover or onset field such 
that it separates the linear regime from the non-linear 
regime along the field axis. Note that temperature does 
not enter explicitly in Eq. ([6]) but does so through the 
temperature-dependent a (T). Thus, according to Eq. 
(J6]) the field scale for nonlinearity is determined solely by 
the linear conductivity <j . This is indeed a very signif- 
icant relation in that the same relation holds good even 
when <t is changed by some other variable like magnetic 
field or simply quenched disorder (see Eq. [5]). Presently, 
a theoretical understanding of this result is lacking but 
diverse disordered systems with localized states^ includ- 
ing amorphous and doped semiconductors, composites at 
low temperatures have been found to obey such scaling, 
all with positive exponents xt- 

At large fields F 3> F a , FV curves often tend to be- 
come independent of temperature. From Eq. ([3]) this fol- 
lows immediately if g{x) ~ x x l XT at large x for xt > 0. 
Thus, at large fields the conductivity varies as a power- 
law with an exponent Zt- 

<j(T, F) ~ F ZT , z T = l/x T for x T > 0. (7) 

This prediction provides a self-consistency check for the 
exponent xt as it may be determined using two inde- 
pendent methods: one from a set of I-V curves at low 
fields and the other from a single I- V curve at high fields. 
This has been also amply verified 4 ^ in systems mentioned 
earlier. It is also clear from this that in case of negative 
Xt one can not expect temperature-independence of con- 
ductivity at large field although the scaling itself can not 
specify any particular functional form. 

Composite systems at room temperature are presently 
the only disordered ones where the nonlinearity exponent 
has been explaine d 47 ' 48 although they are not, strictly 
speaking, comparable to the systems under study in this 
paper. Firstly, conduction in composites (at least) at 
room temperature is not known to take place through 
localized states and is rather metallic since resistance in- 
creases with temperature. Secondly, it is not clear how 
results in a percolating network 49 (which composites are 
supposed to be prototypes of) could be directly applied 
to disordered systems such as CPs since a hopping net- 
work is considered to be a percolating network always at 
criticality 5 - with infinite correlation length. Neverthe- 
less, it may be instructive to review non-Ohmic conduc- 
tion in these systems. In this case, I- V curves have been 



measured at a fixed temperature (i.e. room temperature) 
in samples with various degree of disorder characterized 
by the parameter p, fraction of conductors in a sample. 
Near the percolation threshold p c (i.e. p > p c ), I-V's 
turn out to be nonlinear even at room temperature with 
conductance increasing with the applied bias. As men- 
tioned earlier, the increase in conductance is due to open- 
ing up of new channels of conductions - from tunneling 
across closely spaced tips of pairs of branches of conduc- 
tors dangling off the backbone. The bias scale was found 
to vary as the inverse of the correlation length £ : 

v (p) < r 1 ~ s;/', (8) 

so that the nonlinearity exponent x p is given by x p < v /t 
where v and t are the correlation and conductivity expo- 
nents respectively and E (p) = V = 0) ~ (p — p c ) ■ 
The subscript in the exponent denotes the parameter 
that is varied to change S . v/t is about 0.45 in three 
dimension and consistent with the values of x p obtained 
experimentally in discontinuous gold filrn^ 7 - and carbon- 
wa»^. In a sample of the same carbon-wax system 
with a given p (i.e. fixed disorder), I-V curves were 
measured at different temperatures^ 1 -. The onset volt- 
age V (T) scaled as before with the linear conductance 
S (T) = E(T, V = 0) as 

V (T)~E **. (9) 

Interestingly, xt was found to have the same value as x p 
defined in Eq. ©. 

It is seen from Eqs. © and ([7]) that the nonlinear- 
ity in a system can be characterized by two quantities: 
the field scale F Q that determines the onset and the non- 
linearity exponent xm that determines the degree (i.e. 
steepness of increase of the conductivity with field) . Both 
the scaling function g and the nonlinearity exponent xt 
are obviously determined by the details of the conduction 
mechanism under study. In principle, the problem of non- 
Ohmic conduction is solved if an expression for I(M,V) 
is found where the I- V data are taken at different values 
of the parameter M. Let us now consider several models 
and see how they conform to general scaling formulation. 

2. Variable range hopping under field 

In disordered systems, an increasing electric field aligns 
an increasing number of empty and accessible states to 
the occupied states allowing charge carriers to move via 
phonon- assisted tunneling or hopping transitions. With 
decreasing temperature the mean hopping length Rh 
grows as 

R h = (a/2)(T /T) m , (10) 

where Tq and m are same as in Eq. |T]). With higher 
electric fields the energy eFRh gained by an electron may 
become comparable to ksT and lead to deviation from 
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Ohmic behavior. Theorie o 40 ' 42 ' 43 predict two characteris- 
tic fields Fi and F u such that the non-Ohmic conductivity 
at intermediate fields F < k B T / ea is given by 



a(T,F)=a(T,0) ^p{0) 



(11) 



where er(T, 0) is given by Eq. (p} and L is a length related 
to the hopping length Rh- Clearly, Fi — k B T/eL could 
be identified with the onset field scale F Q discussed above. 
L is proportional to i?^ where /j, is either 1 (Refs. and 
H) or 2 (Ref. S3). Thus, it follows from Eqs. CP and 
(JTUJ) that F a varies with T as a power-law: 



F ^Ft~ T a /aT a -\ a = 1 + m M , 



(12) 



with the same m as in Eq. Q}. According to Eq. (fl2]l, 
the exponent a is always positive and greater than 1. 
In the literature however, authors traditionally focus on 
the temperature dependence of L ~ y-mp rather than 
the field scale. The compliance of this with experimen- 
tal results 5 ^ in amorphous and doped semiconductors is 
rather poor. In many cases the power-law (i.e, a straight 
line in a log-log plot of F a vs. T) is not observed. When 
the plot appears like a straight line, the exponent often 
turns out to be randomly different from the expected 
value in Eq. (fT2"l) albeit with few exceptions. In the 
limit of large fields F > k B T '/ea theorie s 41 ' 42 agree with 
'activationless' hopping at least qualitatively. In such sit- 
uation the energy gained by an electron is large enough 
to hopp without any need to absorb any phonon. As 
a consequence, the conductivity becomes independent of 
temperature and is given by 



cr ~ exp 



F„ 



F„ 



kuT n 



ai- 



(13) 



with the same m and T Q as in Eq. ([T]). a\ is a numerical 
constant equal to unity— when m=l/4. Activationless 
hopping has been observed by many authors 5 — irrespec- 
tive of the value of m. However, inexplicably the partic- 
ular expression of high field conductivity in Eq. (|13[) has 
been found so far only in systems with m=l/2. On the 
contrary, in several systems 5 ^ with m=l/4, the field de- 
pendence is well described, albeit empirically, by a power- 
law a ~ F z reminiscent of Eq. ([T]>. As mentioned earlier, 
Eqs. (IT5|) and @ illustrate the possibility of having self- 
consistent relations involving the same parameters such 
as a, T a . This has led to quantitative disagreemen t 28 ' 54 
in that experimental values of F u turn out to be always 
greater than the calculated ones. 

In presence of coulomb interactions with m=l/2, some 
author a 28 ' 54 ' 55 have suggested that the field- dependent 
conductivity can be written in the following manner: 



a(F, T) w exp 



2R h 



a-l 



eFR h 



AR h T k B T 



(14) 



In the activationless regime the last two terms in the 
exponent cancel each other and Eq. (|13l) is recovered. 



But there is no unanimity in the value of a\ which can be 
1/2^1, l££ or 1.442 8 .. More significantly, it is to be noted 
that while the above expression at least in intermediate 
fields is compatible with the scaling (Eq. [5]), the large 
field conductivity as given by Eq. (fT3j) is not. 



3. Fluctuation-induced tunneling under field 

Kaiser et a\M^- extended the usual FIT model^I at 
small fields by including the backflow current, the tun- 
neling near the top of the barrier and thermal activation 
over the barrier to obtain the following phenomenological 
expression for conductivity: 



ao exp(V/V ) 
1 + h [eMV/Vo) - 1] ' 



(15) 



where <7q is given by Eq. ([3]) and the parameter h = 
co/coo [h < 1), (Too being the value of cr at large voltages 
V. V is a voltage scale factor. At small bias V < V , the 
conductivity increases exponentially with V but at large 
bias V ~S> V , the conductivity increases to the satura- 
tion value Coo. The quantities V and have not been 
related to parameters of the model and their tempera- 
ture dependence remain unknown. The above expression 
is claimed to give a good description of the nonlinear I- 
V characteristics at different temperatures in quasi-lD 
systems like PA nanofibre, single wall carbon nanotube 
network^. However, Yin et al<^ reported systematic de- 
viations between the I-V data and the fitting curves at 
low voltages in several quasi-one dimensional systems like 
individual PPy nanotube, PEDOT nanowire, PA nan- 
otube, and CdS nanopore although overall fits appeared 
reasonable. Moreover, h was found to have a substantial 
dependence upon temperature, increasing by as much as 
two orders of magnitude with temperature in case of PPy 
nanotube. This is contrary to expectations since tJoo is 
expected to be independent of temperature at large bias. 
There have been also several instances^ 8 - where h pa 
indicating that Ooo is probably very large. As discussed 
above, it is the variation of V with temperature or any 
other variable that is of interest here. 



4-. Multi-step tunneling 

Considering the process of multi-step indirect tun- 
neling via n-localized states in a disordered sys- 
tem, Glazman-Matveev (CM) proposed the following 
expression 4 ^ for the conductance through a amorphous 
semiconductor thin film (i.e. tunnel barrier) of thickness 
w under bias V (eV 3> k B T and p n = n — 2/(n + 1)): 



l 



S + EaV 1 - 33 + E 3 V 2 



E 4 V 36 



£ 5 V 



4.67 



EeV" 



',.71 



(16) 
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where Eo = Erf + Si. E^ accounts for the direct tun- 
neling and Ei for the elastic resonant tunneling via one 
localized state. Each term (n > 1) in the series arises 
out of rare events when a number of localized states hap- 
pen to be arranged physically as well as energetically in 
such a way that an electron can traverse a sample length 
via multi-step inelastic tunneling. Each term may be 
thought to constitute a separate channel of conduction 
involving a definite number of localized states. Thus, the 
macroscopic nonlincarity in this GM-model results from 
two contributions - primarily, appearance of a new chan- 
nel with increasing bias and secondarily, nonlinearity of 
each such channel. This multi-step tunneling model has 
been widely invoked to explain relevant data in various 
tunnel junctions^, manganite o 57 ' 58 . Interestingly, the 
manganites samples were not necessarily in the form of 
junctions. Nevertheless, it was thought fit to apply GM- 
model because of possibilities of tunneling across grain 
boundaries, or insulating barriers, separating metallic 
phases within samples. Considering the similar structure 
of conducting polymers i.e. the polymer fibrils consist- 
ing of quasi- metallic lengths with intervening insulating 
regions, CPs are expected to be also candidates for appli- 
cation of the GM-model. It is shown below how a couple 
of assumptions about the coefficients E„ make the GM- 
expression (Eq. I16[) not only compatible with the general 
scaling formulation (Eq. [5]) but also yield the same rela- 
tion between the nonlinearity exponent xt and the bias 
exponent zt as in Eq. ([7]). 

Scaling: While the I-V data at different temperatures 
have been fitted well by Eq. (|16p with few nonlinear 
term o i 58 , no attempt was made to analyze systemati- 
cally the temperature variation of the coefficients E„ in 
Eq. (fl6|) . Note that Eq. (fl6|) as written is not compat- 
ible with the scaling (Eq. [SJ. One way to achieve this 
is to assume that for n > 2, the coefficients satisfy the 
following 



the following relation: 



E„ = c n E Vb p 
so that Eq. p^|) can be written as 



(17) 



E (V 



V 



1.33 



V 

r3{ y 



2.5 



V 



3.6 



C5[ v 



V 

(6[ y 



5.71 



(18) 



The above expression explicitly implies the existence of 
a single bias scale V as in Eq. ([5]). The right hand side 
of the above equation is thus really a scaling function 
as defined in Eq. ([5]). The same holds if conductivity 
and field are used in place of conductance and bias. It is 
seen from Eq. (|17p that the constants c n are determined 
once the coefficients E„ and scale V are known. The ex- 
pression (Eq. [18]) suggests that experimental E-V curves 
at different T's may be collapsed into a single curve by 
suitable choices of E G and bias scale V a at each temper- 
ature. Note that at V = V , the coefficients c„'s satisfy 



E' 



(19) 



If, for example, E(Vo) = 2E Q , we have c « = 1- 

Nonlinearity exponent xt- The GM model does not 
give any guidance on the temperature dependence of the 
coefficients E„ under the condition eV 3> ksT for which 
Eq. pi))) is valid. Therefore, one needs to fall back on 
the empirical relation (Eq. [6]). Alternatively, if the co- 
efficients E„ are known from fittings one can then find 
V in the following manner. Let n = n Q be the lowest 
channel with nonzero E„ o . From Eqs. (jU and (T17J1) the 
onset bias V is given by Y, no V Pn ° ~ E D which leads to 
the expression for V as 



Vo 



_Eo_ 



l/Pno 



(20) 



This is consistent with the more general assumption (Eq. 
rTT|) . In absence of any theoretical guidance, we make 
here a second assumption that E ra follows a power law 
behavior with Eq as given by 



5V" 



(21) 



where y n is an exponent. Incorporating this result in Eq. 
([20| and comparing with Eq. (|6]), one has 



x T 



Pn 



Pn 



(22) 



for any n > 2. The second equality on the right hand 
side of the above equation follows from Eqs. (fT7|) and 
(I2T1) . Notice that while the left hand side of Eq. (|2"2]) 
is a constant, the right hand side carries the index n. 
Thus, if n Q — 2, then xt = |(1 — 2/2)- If n a = 3, then 
xt = I (1 — 2/3) and so on. All three phenomenologi- 
cal relations Eqs. ([6]), (fl"7|). and (f2~Tj) . only two of which 
are independent, underline a basic assumption that non- 
linear scales are determined by the corresponding linear 
conductivity. Eq. (|22[) . a somewhat generalized version 
of the relation considered earlier by Chakrabarty et al4£, 
is remarkable for several conclusions that immediately 
follow from it: 

i) xt is when y rio = 1. It is positive or negative 
depending upon whether y„ o is less or greater than 1. 
Furthermore, its maximum possible value is \/p2 or 0.75 
as p n has the lowest value for n=2. Thus, we have 
-1/2 < x T < 3/4. The lower limit of -1/2 is obtained 
from general arguments^. 

ii) Since the left hand side of Eq. (f22j) is independent 
of n, and p n increases with n, it follows that for xt > 0, 
y n must decrease with n and for xt < 0, y n must increase 
with n. 



1 > V2 > V3 > ■■ > Vn, 
1 < V2 < < ■■ < Vn, 



X T >0 

x T < 



(23) 
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Thus for xt < 0, we have y n — —xtPu for y n 3> 1 at 
large n. If any E„ in the GM expression is zero, the 
corresponding y n is obviously excluded from the above. 

hi) For xt > 0, the above naturally raises question 
about the limit of y n for large n. Considering the fact 
that the conductance of a hopping system at large bias 
may tend to become independent of T or in other words, 
So, the lower limit of y n may be taken as zero i.e. as n — > 
oo, y n — > 0. Note that y n is always a positive quantity 
since the coefficient £„ is predicted to be a decreasing 
function of r*2i. 

iv) Since p n increases with n, Eq. (|22[) requires that 
for xt > 0, the series (Eq. I16[) must terminate at n — 
where the highest channel is given by p nL -i < 1 /xt < 

Pn L - 

Large bias exponent zt- If the field-dependent con- 
ductance of a system continues to be described by Eq. 
((THJ) (e.g. there is no appearance of negative differential 
conductance) it follows that at large bias or field the con- 
ductance varies with bias as a power- law: E rj E„ £ V Pri L 
where til is the index of the highest allowed channel 
(see the point (iv) above). Such a power-law variation 
is also predicted by the scaling consideration (Eq. [7]) for 
xt > with zt = Pn L - Now, from Eq. (|2"2"j) . we obtain 
Xt = (1 - Vn L )/p nL so that zt = (1 - y nL )/x T = l/x T 
following the point (hi). Thus, we obtain zt — 1/^t as 
follows from general scaling consideration (Eq. [7]). 

III. EXPERIMENTAL 

The three systems used in this study were doped PPy 
in powder and film forms, and doped PEDOT in powder 
form. All together five samples (see Table I) have been 
studied. All samples were synthesized using chemical 
polymerization. To aqueous solution of pyrrole monomer 
(0.05 M), 100 ml of aqueous FeCl 3 (0.1 M) solution was 
added dropwise leading to almost instant polymerization. 
Black precipitate of polypyrrole was separated from the 
solution and the resulting powder was then vacuum dried 
and pressed into disk shaped pellets, 8 mm in diame- 
ter and 1 mm in thickness (sample 1). Further details 
can be found in Ref. [24]. A similar procedure was fol- 
lowed for synthesis of doped PEDOT (sample 2) using 
monomer cthylcncdioxythiophene (EDOT) dissolved in 
water with 0.2M dodecylbenzenesulphonic acid. Subse- 
quently, an aqueous solution of the oxidizing agent FeCl3 



TABLE I. Sample growth conditions and parameters (Eq. [TJ- 



System 


Sample 


Oxidant 


<t o (300K) 


m 


T 




No. 




S/cm 




K 


PPy (powder) 


1 


FeCl 3 


5.7xl0~ 4 


1/4 


1.7xlO b 


PEDOT (powder) 


2 


FeCl 3 


2.6xl0" 2 


1/2 


1798 


PPy (film) 


3 


HC1:H 2 02 


10 


1/2 


2340 


PPy (film) 


4 


HC1:H 2 2 


12 


1/2 


950 


PPy (film) 


5 


HC1:H 2 2 


8.5 


1/2 


1600 



-1-0.25/1^-0.25^ 

I (K. ) 
0.1 0.2 0.3 




T 5 (K 5 ) 

FIG. 1. (Color online) Variation of Ohmic conductivity <7o vs. 
T~ m with m = 1/4, 1/2 for four conducting polymer samples 
as shown. For clarity, the PEDOT data have been shifted 
upwards by a factor as shown. Solid lines are linear fits to the 
data. 



was added dropwise on constant stirring under inert at- 
mosphere. EDOT:FeCl3 molar ratio was 1.167:1. Af- 
ter overnight reaction, the precipitate was thoroughly 
washed with ethanol, vacuum dried, and pelletized in 
form of disks of same dimensions as above. PPy films 
(samples 3-5) were synthesized from pyrrole vapor in so- 
lutions of HCl:H2O2=l:1000 while the solvent was cooled 
to 282, 283 and 300 K respectively. Further details can 
be found in Ref. 2_8 (the sample 3 here is the same as the 
sample A in the latter). Thicknesses were 70, 77 and 60 
nm respectively. 

Table I shows typical room temperature conductivities 
of the samples. For transport measurements, thin copper 
wires were attached to both ends of the disks using silver 
paint whereas leads were connected to pre-patterned Pt 
electrodes on the substrate containing PPy filrn^. Two 
and four probe dc transport measurements gave similar 
results indicating negligible contact resistances in sam- 
ples. Low temperature measurements were done in liquid 
helium cryostats in the temperature range 2.1-300 K. The 
PPy and PEDOT samples were placed on sapphire sub- 
strates with Apeizon N-grease and data were taken under 
the constant current condition. In case of PPy films, data 
were taken under the constant voltage condition. Tem- 
perature was stabilized to better than ±50 mK for the 
/- V measurements. Maximum current levels were kept 
low to minimize joule heating in the samples. For exam- 
ple, the maximum current through the PEDOT sample 
at 26 K was limited to 0.22 mA as an increased current 
led to instability due to joule heating. However this be- 
havior was completely reversible and sample resistance 
returned to the initial value after removal of current. 
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IV. RESULTS 

Fig. 1 shows Ohmic conductivities a vs. T~ m for 
four samples - one each of PPy (powder) and PEDOT, 
two of PPy (film) . Similar data for sample 3 have been 
already presented in Ref. [28]. Excellent linearity in the 
data shows that all the samples obey variable range hop- 
ping conduction in the range of temperatures measured. 
While m is 1/4 in PPy (powder), it is 1/2 in all sam- 
ples of other two systems. There is a transition at 30 K 
in PPy films from ES-VRH mechanism to activated pro- 
cess. A similar transition from ES-VRH to M-VRH also 
takes place in PEDOT at about 38 K. All the samples 
were in insulating regimes. To's obtained from slopes of 
the fitted lines are of the order of 10 3 K and are given 
in Table I. The exponent 1/2 has been shown to result 
from coulomb gap at the fermi level 3 - . Parameters (i.e. 
m, T ) in PPy (powder) and PEDOT agree with those 
reported earlie r 24 ! 25 . Values of selected fj (T)'s are given 
in Table II. It is seen from Table I that of the first three 
samples, the sample 1 (PPy powder) has the least con- 
ductivity at room temperature while the sample 3 (PPy 
film) has the highest. We discuss the systems in this or- 
der as the maximum normalized conductivity (see below) 
also follows the same order. 



A. Polypyrrole (powder) 

The field dependence of dc conductivity a of a 
PPy(powder) sample at different temperatures ranging 
from 80K to 300K are shown in Fig. 2. The non-linear 
response of conductivity to the application of electric field 
can be clearly seen from the figure. The sample had a 
zero-bias linear conductance So of 0.003 S at room tem- 
perature. A typical behavior at a constant temperature 
is that the conductivity remains constant for small fields 
and then starts increasing with increase in field. The 
value of the field at which conductivity starts deviating 
from its linear value a is the onset field F Q . A criterion 
to determine the latter is discussed below. With fur- 
ther increase in the bias, the conductivity continues to 
increase monotonically. As temperature is decreased, a a 
decreases and the sample seems to become nonlinear at a 
field greater than the one required at a higher tempera- 
ture i.e. F Q increases with decreasing temperature. This 
behavior is opposite to those found in other two systems 
(Figs. 4 and 6). Data at each temperature were first 
fitted to Eq. (fl"5|) with all three parameters <r , V and h 
being free but fits were rather poor. In contrast, fits to 
the GM-expression (Eq. II 6|) containing terms up to n=5 
were reasonably good upto cr/a ~ 2.5 and are shown 
in Fig. 2 (lines). But as seen particularly at T=150 
and 195 K, conductivity increased faster than the fitted 
curves at higher bias and could not be accounted even by 
including a n=7 term. This may be related to incipient 
negative differential conductance regime as is evident in 
similar data in Fig. 2c of Ref. l26l . Surprisingly, nonlinear 



1 o 3 _ PPy (powder) 
: 300K ^hm- 





' ' ' ' 

10 3 10 2 10 1 10° 10 1 
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FIG. 2. (Color online) Variation of conductivity vs. elec- 
tric field in a doped Polypyrrole pellet (sample 1) at different 
temperatures as indicated. The dotted line schematically in- 
dicates the movement of the onset bias with increasing tem- 
perature. The solid lines are fits to the Glazman-Matveev 
(GM) expression (Eq. I16[) . See text for details. 



least square fittings at all temperatures led to either very 
small or negative values for coefficients of even terms (i.e. 
n=2,4) so that final fittings were done using only three 
terms: a(V) =a Q + a 3 V 5 / 2 + fj 5 V 14 / 3 . Fitted values of 
the coefficients are given in Table II. 

Fig. 3 shows the result of making the data in Fig. 2 
collapse into a single curve by suitable scaling. It is con- 
venient to start with a temperature such that the data at 
that temperature are predominantly Ohmic but contain 
minimum non-Ohmic regime. In the present case the ap- 
propriate starting temperature is 80 K. The conductivity 
was scaled by its Ohmic values a . For the field, any 
arbitrary choice (e.g. 1) for F Q would do as far as data 
collapse is concerned. For the next higher temperature, 
the conductivity was scaled as before but F Q was adjusted 
in such a way that this set of data merged with the earlier 
one as best as possible. The same procedure was then re- 
peated for all the other temperatures in increasing order. 
Note that in this method F Q is determined only upto a 
constant value. Multiplying all F 's by a constant only 
shifts the merged curve along the field axis without alter- 
ing the curve anyway. To facilitate comparison, and for 
F a to be interpreted as an onset field, its scale was fixed 
by adopting an uniform criterion that the conductivity 
at the onset field would be double of its ohmic value i.e. 
cr(F ) = 2rr(0) = 2cr . The excellent data collapse up to 
about a / a « 3 seen in Fig. 3 proves the existence of a 
field scale at each temperature. F a thus obtained follow- 
ing the above criterion is plotted with log-log axes as a 
function of both temperature T (closed symbol) and the 
corresponding a Q (open symbol) in the inset. The solid 
line indicates a power law F a ~ er|; T with an exponent 
Xt being equal to —0.329 ± 0.014 which is negative as 
suggested by the orientation of the dotted line in Fig. 1. 
However, no such relation is apparent in the functional 
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FIG. 3. (Color online) Scaling of the same data as in Fig. [2]to 
achieve data collapse as shown. The scale of the onset field F 
is fixed by adopting a common criterion (j{F ) = 2a in this 
paper. The solid line is a fit to the scaled version (Eq. I18p of 
GM expression with only two odd (n=3,5) nonlinear terms: 
a I (Jo = I + 0.54g 5/2 + 0.42g 14/3 (g = F/F a ). The dashed line 
is a plot of the same expression without the n=5 term. Inset 
shows two log-log plots of F vs. <r (open symbol) and T 
(closed symbol) . The solid line is a linear fit to the data with 
a slope xt as shown. 



dependence on T. Results thus validate the scaling as 
given in Eqs. (O and ©. Since all the curves in Fig. 2, 
that were fitted by the GM-expression also collapse on to 
a single curve it is expected that the latter would also be 
fitted by the scaled form (Eq. [T5)) . This is indeed con- 
firmed by the solid line in Fig. 3, which is a fit according 
to a/a Q = l + c 3 (F/F ) 5/2 + c 5 (F/F ) u/3 with c 3 = 0.54 
and C5 = 0.42. The relative high value of C5 confirms 
the rapid increase of the conductivity with field. A plot 
(dashed line) of the same fitting expression without the 
n=5 term is also shown to highlight the contribution of 
the higher nonlinear term. The divergence of the data 
from the fitted curves at higher fields (<j/<j ~ 3) as seen 
in the figure is due to th negative differential conductance 
as discussed above. c„'s in the fitting expression are sim- 
ply averages of c„'s calculated at each temperature from 
Eq. (fT7| using cr„'s and V — 0.1F o and given in Table 
II (for calculation of C5, data at 80 and 100 K were ig- 
nored as non-Ohmic regimes were small). The highest 
nonlinear term in the fitting expression depends upon 
the maximum value of the measured normalized conduc- 
tivity, (cr/(Jo)max which is about 4 in the present case. 
Notice that C3 + C5 = 0.96 that is close to 1 as expected 
from Eq. (|T9|) by applying the criterion for F a which also 
leads to a condition for the scaling function in Eq. ([S]): 
g(l)=2 seen also in Fig. 3. The deviation of the sum 



of c„'s from 1 is an indicator of how well the criterion 
for the field scale was implemented during scaling and 
quality of overall scaling. 

It may be mentioned here that if one is interested only 
in the nonlinearity exponent and not the scaling func- 
tion, it could be obtained from experimental data by an- 
other method used by Gefen et al. 4 - for characterizing 
the crossover to the nonlinear regime in a percolating 
system. In this method, one defines the crossover field 
F Q such that the conductance at this field deviates by 
an arbitrarily chosen factor e from its zero-field value i.e. 
S(_F ) = S (l + e). Obviously, the value of the exponent 
should not depend upon the choice of e as verified in dis- 
continuous gold films^. Two values of e, 0.1 and 0.4, 
were considered but no significant variation in the expo- 
nent was observed, e = 1 coincides with the criterion 
adopted in this paper. Clearly, this method will work 
as long as S Q can be obtained from data. But particu- 
larly at low temperatures, Ohmic conductivities become 
too small to be above the noise floor of measurements. 
Hence, measurements are feasible only in highly nonlin- 
ear regimes and consequently, S 's can not be obtained 
directly from the data. In such cases, the method of scal- 
ing provides an alternative way to take into account such 
nonlinear data as illustrated in the cases discussed next. 



B. PEDOT (powder) 

The field dependence of dc conductivity a of a PE- 
DOT(powder) system at different temperatures ranging 
from 4.9 to 26 K are shown in Fig. 4. The nonlinear re- 
sponse of conductivity to the application of electric field 
can be clearly seen from the figure. The sample had an 
Ohmic conductance T. of 0.026 S at room temperature. 
It is observed from the figure that the dotted line points 
towards right in contrast to the behavior in the previ- 
ous case i.e. the onset field F Q increases with increas- 
ing temperature or conductivity. In fact, the orientation 
of the dotted line determines the overall shape of the 
curves which are somewhat divergent with bias in PPy 
(powder) (Fig. 2) but in contrast, appear convergent 
in PEDOT (powder) (Fig. 4). The converging feature 
naturally indicates eventual temperature independence 
of conductivity at large fields. As before, data at each 
temperature were fitted to the GM-expression (Eq. \W\i 
containing terms up to n=5. In this case, the least square 
method at all temperatures led to either very small or 
negetive values for coefficients of odd terms (i.e. n=3,5) 
so that final fittings were done using only three terms: 
a{V) = cr Q + a 2 V 4 ^ 3 + a 4 V 18 / 5 . The Ohmic conductivi- 
ties at 4.9 and 6.5 K were obtained from the extrapolated 
line, In a a vs. T~°' 5 in Fig. 1. The bias ranges at 10 and 
12 K were too limited to yield reliable values for C4 S. Fits 
(lines in Fig. 4) are seen to be excellent. Fitted values of 
the coefficients are given in Table II. 

All the curves in the figure, starting with the one at 
26 K and using the same criterion for F a as in the pre- 
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FIG. 4. (Color online) Variation of conductivity vs. electric 
field in doped PEDOT (sample 2) at different temperatures as 
indicated. The dotted line schematically indicates the move- 
ment of the onset bias with increasing temperature. The solid 
lines are fits to the GM-expression (Eq. I16|) . See text for de- 
tails. 



T(K) 

1 10 





■ 3 


I X T= 


0.158±0.003 




PEDOT(powder) 


100 


r2 




/y 




I OV/cm) 






' 1 






if/ 

LL ! 1 I 


c 

.b 

10 




' J 


j j ■■■ 


..J ■ 




: 10" 10 


10" 8 10" 6 
a (S/cm) 


10 " 4 IJ/++ 4.9K: 

If > 6.5 ; 












f O 10 . 












jf O 12 . 












v 16.5 


1 






ivo ipvQOI»»^ 
i i 




O 26 - 
■■i i ...... 



10" 3 10" 2 10" 1 10° 10 1 
F/F o 

FIG. 5. (Color online) Scaling of the same data as in Fig. [4] 
to achieve data collapse as shown. The solid line is a fit to 
Eq. (|18[) with only even (n—2,4) nonlinear terms: a/a = 
1 + 0.92g 4/3 + 0.011<2 18/5 (q = F/F a ). The dashed line is a 
plot of the same expression without the n=4 term. The dash- 
dotted line is a fit to an exponential function (Eq. Hip . Inset 
shows two log-log plots of the scaling field, F vs. a (open 
symbol) and T (closed symbol). The solid line is a linear fit 
to the data (open symbols) with slope as shown. 



vious case, could be made to collapse into a single curve 
as shown in Fig. 5. Since the data for T < 7 K have 
no linear region, it may be noted that the method in 
Ref. [47| for determining the exponent can not be ap- 
plied. However for scaling purpose, extrapolated <7 's 
were used. However, when conductivities are scaled with 



these rj 's, the scaled data were all larger than the rest 
of data at higher temperatures. Absence of any over- 
lapping data sets may leave some margin of error in F a 
which needs to be adjusted to complete the process of 
data collapse. The larger the gap between data sets is, 
the higher the margin of error is. The margin of er- 
ror may be considerably reduced when, as in the present 
case, the collapsed curve is compared in somewhat self- 
consistent manner to some function it is expected to fol- 
low. The solid line in the figure is an excellent fit to 
a [<x = l+c 2 (F/F ) i/3 +c 4 (F/F ) 18/5 with c 2 = 0.99 and 
C4 = 0.015 so that c 2 +C4 = 1.005 as expected. c 2 is the 
average of c 2 's (with data at T=26 K ignored) whereas 
C4 is consistent with the values of limited number of C4's. 
A plot (dashed line) with only n=2 term is also shown to 
highlight the contribution of the higher nonlinear term, 
which, as seen, is quite significant in this case although 
C4 is quite small compared to c 2 . Note that (a/a )max 
is about 100 compared to 4 in the previous case. In- 
set shows log-log plots of F a vs. temperature T (closed 
symbol) and the corresponding a Q (open symbol). The 
solid line through open symbols indicates a power law 
with an exponent of 0.158 ± 0.003. The exponent has a 
positive value in accordance with the orientation of the 
dotted line in Fig. 4. No reasonable straight line could 
be drawn through closed symbols. 



C. Polypyrrole (film) 

Two sets of field-dependent conductivities of PPy films 
are presented in Fig. 6. The panel a shows data taken 
in the sample 3 at different temperatures as marked and 
is similar to Figs. 2 and 4 whereas the panel b shows 
data at Tm 20 K taken in three different samples as in- 
dicated. In the former, the initial conductance changed 
due to change in the temperature whereas in the latter, 
the same was achieved by having different quenched dis- 
order in the samples. As seen in the Figs. 2 and 6, the 
basic qualitative response to the electric field is similar in 
all the PPy samples irrespective of structure (i.e. powder 
or film) or sample condition in that at any temperature 
the conductivity behaves as a monotonically increasing 
function of field, starting from a constant value at small 
bias. However, a closer look reveals subtle differences 
as illustrated by the orientations of the schematic lines 
which indicate the movement of the onset field F a with 
increasing linear conductivity. In Fig. 6a, the line points 
towards right (i.e. F a increases with the linear conduc- 
tivity) as in Fig. 4 whereas in Fig. 6b, it points towards 
left (i.e. F Q decreases with the linear conductivity) as 
in powder (Fig. 2). Note that although both panels 
have one common set of data (sample 3 at T=20 K) this 
did not prevent two sets of data with two different driv- 
ing variables - namely, temperature and disorder - from 
exhibiting opposite behavior in the onset field F a . As 
mentioned in the previous case, the relation of the over- 
all shape of the curves with the orientation of the dotted 
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FIG. 6. (Color online) (a) Variation of conductivity vs. field 
in a doped PPy film (sample 3) at different temperatures as 
indicated. The dotted line schematically indicates the move- 
ment of the onset field with increasing temperature. The solid 
lines are fits to the GM expression (Eq. I16[) with parameters 
given in Table II. (b) Similar data as on the left panel taken 
at the same temperature 20 K in three different samples as 
indicated with different quenched disorder. The dotted line 
schematically indicates the movement of the onset bias with 
increasing conductivity. The solid lines are fits using the GM 
model. See text for details. Note that data represented by 
solid diamonds in both panels are same. 



lines is now clearly seen in the two panels - the left one 
showing a convergent behavior as in Fig. 4 and the right 
one a divergent behavior as in Fig. 2 - in the same sys- 
tem. At the low temperatures (2-5 K) conductivity in- 
deed become nearly independent of temperature at large 
biases as seen in the panel a - a feature that was strongly 
hinted in the data of PEDOT. The same is also appar- 
ent from roughly same values of o§ at low temperatures 
(see Table II). As before, data at each temperature or 
disorder were fitted to the GM-expression (Eq. [T5|) con- 
taining terms up to n=6. In this case, the least square 
method at all temperatures led to either very small or 
negative values for coefficients of odd terms (i.e., n=3,5) 
so that final fittings were done using only four terms: 
a(V) = a + a 2 V 4 / 3 + a 4 V 18 / 5 + <7 6 V 40 / 7 . Fittings to 
data at T < 7 K used the Ohmic conductivities obtained 
from extrapolation of In a Q vs. T~ 5 line. 

Fig. 7 shows scaling of the data in two panels of 
Fig. 6. There are two curves showing data collapse - 
one (labeled b) belongs to different samples in panel b 
and another (labeled a) belongs to the sample 3 at dif- 
ferent temperatures in panel a. Data collapse has been 
achieved by following the same procedure as adopted in 
previous cases including the criterion for fixing the on- 
set field F , namely cr(F ) = 2a . Excellent data col- 
lapses are seen to have been achieved in both cases, with 
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FIG. 7. (Color online) Scaled conductivity a/cro vs. scaled 
field FjF of doped PPy films in two panels of Fig. 6. The 
scaled data b (lines) belonging to different samples are shifted 
upwards for clarity. The solid line is a fit to Eq. (|18[) with 
only two even (n=2,4) nonlinear terms: a ja — l + 0.97g 4 ^ 3 + 
0.0013(? 18/5 (q = F/F ). The scaled data a with symbols be- 
long to different temperatures. The solid line through it is a fit 
to the GM expression (Eq. I18[) with three even (n=2,4,6) non- 
linear terms: a/a a = 1 + 0.97g 4/3 + 0.0024g 18/5 + 6.5 - V° /7 . 
The dashed line is a plot of the same expression but with 
terms up to n=2 only. Similarly, the dotted line is a plot 
with terms up to n=4 only. Inset shows three log-log plots 
of F vs. g (open symbol and cross) and T (closed symbol). 
The crossed symbols, belonging to three different samples, 
have been appropriately shifted (xlO -5 , x0.2) to be within 
the scales as shown. The solid lines are linear fits to the data 
with slopes xt (open symbol) and xd (cross) as shown. 



curve a covering nine orders of magnitude in conductiv- 
ity and five orders in field. The data of curve b have 
been shifted upwards for clarity. Moreover lines, instead 
of symbols, are used in this case to highlight the quality 
of data collapse. In case of curve a, extrapolated <t 's 
were used for scaling particularly at low temperatures. 
The solid lines in the figure are excellent fits to a ja D — 

l + c 2 (F/F f 3 +MF/F ) 1S/5 + c 6 (F/F o ) 40/7 . The co- 
efficients C2, C4 and ~c% are 0.98, 0.0024 and 6.5 x 10~ 8 for 
curve a, and 0.97, 0.0013 and for curve b. In both cases, 
C2 + C4 + CQ k 1 as expected from Eq. (fl9|) . c 2 's are same 
as the average of C2's in Table II (with data at T=2.1 
K ignored) whereas c 4 's are consistent with the values 
of limited number of available C4S. As with PEDOT, 
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the values of cVs are much greater than those of other 
coefficients although the latter's contributions are signif- 
icant as indicated by the dashed and dotted lines. The 
latter are plots with up to n=2 (dashed) and 4 (dotted) 
terms only. This is, of course, partly due to the crite- 
rion adopted for F Q . Note that (a/a ) ma x is about 10 9 
(curve a) compared to 100 and 4 in the previous cases. 
Inset shows log-log plots of F Q thus obtained from scaling 
as function of both temperature T (closed symbol, curve 
a) and the corresponding a (open symbol, curve a and 
cross, curve b). The solid line verifies the power law (Eq. 
[5]) with the exponent xt — 0.228 ± 0.008. No reasonable 
straight line could be drawn through the plot of F D vs. T 
(closed symbols). In case of scaling with disorder, there 
are only three points that yield a tentative value of the 
exponent xu ~ —0.30 (the subscript stands for disorder). 
As expected, xjj compared to xt, has a negative sign. 



D. Other systems in literatures 

Nonlinearity exponents xt in the three systems dis- 
cussed above are shown in Table III. The latter also con- 
tains, for sake of comparison, results from digitized data 
of three other CP systems available in literatures. All 
these systems invariably exhibit the property of scaling 
(Eq. [SJ. The systems include a p-Toluensulfonate (PTS)- 
doped PPy (PPy(R)) nlm 2 ^, PTS-doped polydiacetylene 
(PDA) single crystal^! and iodine-doped polyacetylene 
(PA) nanofibre^l The PPy(R) him (0.1-0.15 mm thick- 
ness) was obtained by electro-deposition at the current 
density of 0.2 mA/cm -2 . The scaled curve of conductiv- 
ities at four temperatures (16-31 K) in this him looked 
very similar to that in Fig. 3. It also yielded a nega- 
tive exponent, -0.16 compared to -0.33 obtained in PPy 
(powder). The temperature- variation of the Ohmic con- 
ductivity of this system has been mentioned in the intro- 
duction. 

PTS-doped PDA crystals are quasi-lD in nature, con- 
sisting of weakly coupled linear parallel chains of co- 
valently bonded carbon atoms. It follow a VRH-type 
conduction (m=0.65-0.70, T D = 2570 K) at low tem- 
peratures with a crossover at about 50 K to activated 
conduction at higher temperatures. Five temperatures 
between 50 and 14 K were used such that correspond- 
ing data could be digitized with some reasonable accu- 
racy from linear current- held plots (viz. Fig. 7 of 
which, particularly at low bias, are more prone to dig- 
itizing errors than conductivity-held ones. Nonetheless, 
Fig. 8 shows clear evidence of scaling in conductivities 
at different temperatures covering nearly six decades in 
conductivity and four decades in field. The solid line 
in the figure is an excellent fit over the whole range to 
ct/ct = l + c 2 (F/F ) 4/3 +c 4 (F/F ) 18/5 . The coefficients 
c 2 and c 4 were 0.98 and 0.000123 respectively and add 
up to nearly 1 as expected. These values were consis- 
tent with the coefficients obtained from fittings of data 
at each temperature to the GM-expression (Eq. Ho]) . The 
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FIG. 8. (Color online) Scaling of the data in a PDA single 
crystal from Ref. [23 to achieve data collapse as shown. The 
solid line is a fit to Eq. (|18[) with only even (n=2,4) nonlinear 
terms: a/a = 1 + 0.98<? 4/3 + 0.000123<? 18/5 (q = F/F ). The 
dashed line is a plot of the same expression without the n=4 
term. Inset shows two log-log plots of the scaling field, F vs. 
a (open symbol) and T (closed symbol). The solid line is a 
linear fit to the data (open symbols) with slope as shown. 



dashed curve is a plot with up to n—2 term only. It was 
found that fitted values of <7 's (also used in scaling) were 
progressively less than those given in the paper with de- 
creasing temperature, by as much as an order of magni- 
tude at 20 K. (a/cr ) max is about 10 5 compared to 10 9 
in PPy film (this work). The nonlinearity exponent was 
xt = 0.51 ± 0.02. It is seen in Fig. 8 that curves partic- 
ularly at lower temperatures tend to rise less rapidly at 
higher fields than lower ones. This is because those por- 
tions of the curves follow F -1 / 2 dependence (Eq. [TBI and 
are outside of the scaling domain^. In fact, the collapsed 
curve really represents an envelope of all scaled curves at 
different temperatures. Such F~ x / 2 dependence in PPy 
film of this work was insignificant (see Fig. 7) although 
both these systems had the VRH exponent m equal to 
1/2- ^ 

Individual iodine-doped PA nanofibrcs of diameters 10- 
40 nm were investigated for transport properties. At low 
temperatures below 30 K, the I- V characteristics follow 
the Zener-type tunneling, £ = S D exp(— F u /F) where F u 
in the exponent depends on the magnitude of the energy 
gap and the effective mass of tunneling electrons^. The 
I- V curves were non-Ohmic and temperature indepen- 
dent up to 30K above which I- V curves were temperature 
dependent. Upon comparing with Eq. (fl~3"|) it becomes 
evident that the above equation really characterizes the 
high-field regime that becomes apparent particularly at 
low temperatures (see also Fig. 6) and is not compat- 
ible with the scaling (Eq. [S]). At higher temperatures 
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TABLE II. Parameters in the GM-expression (Eq. I16[) fitted to conductivity data at various temperature in five samples of 
three CP systems. a n (n > 2) is in unit of S/cm y n_ 2 /(n+ 1 ) < Cn ' s are constants defined by Eq. (|17|l . A blank space indicates 
that the term was not used in the fit. A hyphen means that the data range did not permit reliable determination of the term. 
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the Ohmic conductance exhibits an activated-type con- 
duction and /- V curves again exhibit the same scaling 
behavior shown in Fig. 9 as in other cases discussed 
here. I-V data were obtained in a single iodine-doped 
PA nanofibre of diameter 20 nm at various tempera- 
tures (Fig. 1 of j3_4]). Out of those, data at five temper- 
atures (234-103 K) were found suitable for digitization 
and presented in scaled form in Fig. 9. In spite of some 
inherent noise in the data, high quality of data collapse 
is quite evident. The solid line in the figure represents 
an excellent fit to a/a a = 1 + (F/F ) 4 ^ 3 thus requiring 
only n=2 term in the GM-expression. Authors^ fitted 
the I-V curves to the expression (Eq. IT5j) . The dashed 
curve is a plot of the Kaiser expression (Eq. [15]) with 
/j=0.21. The particular value of h was chosen to satisfy 
the criterion for the bias scale, S(V^) = 2E Q . The fit 
runs above the data for V/V < 1, and runs below the 
data for V/V > 1. Moreover, the variation in h with 
temperature and the saturation at large bias clearly re- 
veals inadequacies of the Kaiser expression (Eq. 1151) in 
describing the field-dependent conductance. Discrepan- 
cies at low bias were already noted^. Inset shows two 
plots of V vs. S with one set of V 's (closed symbols) 
obtained from scaling and another set of V 's (open sym- 
bols) obtained from fittings to Eq. (|15[) . Comparison of 
the two sets shows further limitations of the Kaiser ex- 
pression. The nonlinearity exponent was determined to 
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FIG. 9. (Color online) Scaling of the data in a PA nanofibre 
from Ref. [34|to achieve data collapse as shown. The solid line 
is a fit to the GM-expression (Eq. Ql]): a/a a = 1 + (V/Vb) 4/3 
. The dashed line is a plot of the expression (Eq. I15|l with 
/i=0.21. Inset shows two log-log plots of the scaling bias V vs. 
E . One set of data (open symbols) is obtained from fittings 
to Eq. (|15[) and another set (closed symbol) from scaling. 
The solid line is a linear fit to the data (closed symbols) with 
slope as shown. 
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be x T = 0.33 ±0.02. 



V. DISCUSSIONS 

The last section dealt with four different CP systems, 
namely PPy, PEDOT, PDA and PA in as many forms, 
namely powder, film, crystal and nanofibre. Even in such 
diverse conditions, the scaling phenomena as embodied 
in Eq. ([3]) and demonstrated in Figs. 3,5,7,8 and 9 
stands validated in clear and unambiguous manner. Fur- 
thermore, the figures also confirm the remarkable fact 
that there exists a single field scale in any given sam- 
ple at least within the experimental ranges of field and 
conductivity spanning more than five and nine orders of 
magnitude respectively. This is contrary to the predic- 
tions in the field-dependent VRH theories discussed in 
section II. Thus, an important objective of this paper, 
as stated in the introductory section, is fulfilled. We be- 
lieve that scaling phenomena observed in CP's here, and 
in amorphous- and doped-semiconductors^ earlier, thus 
indicate a general and fundamental property of the class 
of disordered systems with localized states. The scaling 
analysis here follows that of critical phenomena in ther- 
modynamic phase transitions^- and hence, the method 
of analysis described in section III naturally differs from 
that hitherto adopted e.g. in Refs. [H, H3 and HH The 
method, in absence of a proper theory, is primarily phe- 
nomenological but finally yields a concrete number - the 
nonlinear exponent - as a characterization of the underly- 
ing conduction mechanism. Let us now consider details of 
scaling, namely the scaling variables (which F D depends 
on) and scaling function, g and their possible connections 
to the microscopic picture. 



A. Field scale F and scaling variable 

To start with, let us note that at low fields, the vari- 
able range hopping as the conduction mechanism is often 
observed in conducting polymer systems (barring sys- 
tems like PA nanofibre). VRH is actually a phonon- 
activated hopping between localized states irrespective 
of the presence or absence of polarons. In such situa- 
tions, the conductivity as a function of temperature is 
generally described by Eq. ([T]) with m ranging from 1/4 
to 1. Discussed in section II, traditional theories incor- 
porating field effects in VRH conduction have a number 
of predictions or implications which are at variance with 
experimental results: 

First, two field scales F a and F u (corresponding to the 
two length scales, the hopping length Rh and localiza- 
tion radius a respectively) are predicted whereas only 
one scale is experimentally observed; Second, F a is basi- 
cally set by the temperature scale and supposed to vary 
as F D ~ T a (Eq. [T2")) where a = 1 + m/x is a positive 
number and always greater than 1. But log- log plots of 
F Q vs. T (insets in Figs. 3,5,7 and 8) generally devi- 



ate from linearity in varying degrees. They seem to be 
better described by two power-laws with two exponents, 
a\ at low temperatures and cvi at higher temperatures 
with «i > 0i2- Nonetheless, if the data are still subjected 
to linear fittings, slopes (i.e. a's) turned out to be 0.98, 
2.24 and 4.1 60 in PEDOT, PPy (film) and PDA respec- 
tively compared to predicted values of 1.5(2), 1.5(2) and 
1.68(2.36) for fi=l(2) respectively; Third, the problem of 
temperature directly determining the field scale is rather 
dramatically highlighted by the measurements in three 
samples of PPy film (Fig. 6b) at a temperature of 20 K, 
each one having a different quenched microscopic disor- 
der and characterized by the Ohmic conductivity a . If 
F D is indeed set by temperature alone, it should be basi- 
cally same for each sample in the figure as measurements 
were performed at the same temperature. However, In- 
set in Fig. 7 clearly shows that F Q varies with a and 
within the limited range, is compatible with Eq. © al- 
beit with a different nonlinearity exponent. It proves that 
a field scale is an intrinsic property of a disordered sam- 
ple and is a function of various parameters like disorder, 
temperature etc.. Significance of the role of disorder be- 
comes further apparent when compared with the results 
obtained by Bufon et ali^S from measurements at differ- 
ent magnetic fields. Fig. 11 displays a- V curves (dashed) 
at three magnetic fields B as indicated. The conductiv- 
ity generally decreases with increasing magnetic field. It 
is found that the three curves could be made to collapse 
into a single curve (solid) by simply scaling the conduc- 
tivity at each magnetic field by a factor A (arbitrary up 
to a constant factor) as shown. The bias did not need any 
scaling. This means that the field scale was independent 
of B unlike disorder, and consequently, the corresponding 
nonlinearity exponent xb is 0; Fourth, in some systems 
like PPy (powder) (inset of Fig. 3) and PPy(R) film (Ta- 
ble III), a (as well as the nonlinearity exponent xt) is 
negative. This is a serious discrepancy as it is irreconcil- 
able with the theories. 

In contrast, log- log plots of F a vs. a a are consistently 
linear as seen in Figs. 9 and 10, the latter displaying data 
from different CP systems (Figs. 3,5,7 and 8), including 
PPy(R) in the same scales (PA nanofibre is excluded be- 
cause of variables with different units). Linearities in the 
plots give credence to the empirical power-law F Q ~ <j XT 
(Eq. which is radically different from Eq. (fTTj) in that 
the variable in the power-law is not the temperature it- 
self but the linear conductivity a which is, of course, a 
function of temperature. All these considerations suggest 
that cr should be considered as the appropriate scaling 
variable in Eq. ([5]) which then reads as 

cr(M,F)/a = g(Fa ~ XM ), (24) 

where M stands for the variable(s) used to vary a = 
a(M,F = 0). Eq. (|2"3|) allows description of the scaling 
of the field-dependent conduction along various paths in 
the variable space in a natural fashion without any obvi- 
ous contradiction. In analogy with the scaling formula- 
tion of thermodynamic critical phenomena^,, a = de- 
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FIG. 10. (Color online) Field scale vs. linear conductivity in 
various systems. The plot for a-Ge is from Ref. |45| . The data 
of PPy(R) film have been shifted to the left by a factor of 10 4 
to lie within the scales. 
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FIG. 11. (Color online) Conductivity vs. bias voltage (dashed 
curves) in the PPy film (sample 3) at three magnetic fields 
(after Ref. l28l ). Only conductivity of each curve need to be 
scaled with a factor A to collapse on to the solid curve. The 
non-linearity exponent xb is zero. 



fines a 'critical' point. The Ohmic conductivity a Q plays 
the traditional role of temperature T in that it sets the 
field scale which, in turn, must be given by some phys- 
ical length scale. The latter is yet to be explained but 
must be distinct from either a or R^. Similar problem 
of unsolved length scaled exists also on the metallic side 
of the metal-insulator transition 23 in doped semiconduc- 
tors (e.g. Si:P). The linear conductivity of a sample at 
K on the metallic side goes to zero at the transition 
as a power-law as the doping level is decreased to a crit- 
ical value. Such critical behavior is associated with a 
diverging length scale that still remains unsolved. How- 
ever, one must be careful not to carry the analogy with 
temperature too far as the linear conductivity is a non- 
equilibrium quantity which may behave differently than 
the equilibrium quantities 62 . It may be noted that the 
right hand side of Eq. © could still contain a slow vary- 
ing function of temperature or some other variable but 
the dominant variation would still be given by the power- 
law. The doping level like the conducting fraction p in 
composites is the natural quantity for parameterization 
of quenched disorder in CP's. However, in practice the 
doping level often may not be known. Under this situa- 
tion, we use resistivity p (T) = 1/<j (T) as a measure of 
disorder. This is intuitively satisfying but has the draw- 
back of being temperature-dependent. 

While allowing a view of data ranges in six samples 
at a glance, Fig. 10 also provides a basis for quantita- 
tive comparison of the field scales across those systems 
as well as amorphous germanium (a-Ge)^ included for 
reference. The comparison, of course, is not straightfor- 
ward as the slopes xt have different values. However, it 
is seen that CP systems are generally electrically 'soft' 
compared to a-Ge since the onset field in a CP system 
is smaller than that in a-Ge in the displayed range of 
conductivity. For a quantitative comparison, let us con- 



sider a-Ge and PEDOT since both systems have nearly 
same slope (nonlinearity exponent) of 0.16. The prefac- 
tor At in Eq. ([6]) which represents the strength of a field 
scale is 331 and 0.2 (kV/cm)(S/cm)~ 016 respectively i.e., 
the onset field in a-Ge is roughly 1600 times that in PE- 
DOT. Now, simple dimensional considerations lead to an 
expanded expression for F Q : 

Fo = Cl ^(^) XT , (25) 

where C\ is a dimensionless unknown constant, e 2 jh is 
the quantum of conductance. Accordingly, At ~ T Q /a. 
Taking localization lengths in the two systems to be of 
same order of magnitudes (~1 nm) and with T ss 10 8 
and 1800 in a-Ge and PEDOT respectively, the ratio 
F (a-Ge)/F (PEDOT) is about ~ 10 5 , two orders of 
magnitude greater than 1600. VRH theories fare even 
worse. According to Eq. ([I2"]t. F Q ~ l/aT " 1 so that the 
ratio F (a-Ge)/F (PEDOT) is about lSOO^/lO 8 / 4 ~ 1, 
three orders of magnitude less than the experimental 
value. 



B. Scaling function 

At least three forms of nonlinear /- V curves in hop- 
ping systems have been discussed in section II. The field- 
dependent VRH theorie s 40 ' 42 ! 43 predict a simple expo- 
nential (Eq. [TT]) at low and intermediate fields, to which 
the scaled data of PEDOT(powder) have been fitted for 
illustration purpose (dash-dotted line in Fig. 5). The fit- 
ted curve appears to match data well at least numerically 
in the low field region (F/F a ~ 1) but gradually deviates 
(rising faster than data) from it in the high field region as 
predicted. Nevertheless, as shown in the inset of Fig. 5, 
F Q fails to follow the temperature-dependence as given by 
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Eq. (IT21) . Such inconsistent behavior of the VRH theory 
is quite typical. Inadequacies of the Kaiser expression 
(Eq. [T5j) in case of PA nanohber have been discussed 
earlier in details (see Fig. 9). Interestingly, fits^i to Eq. 
(|15p were visually quite acceptable. This brings us to the 
GM-expression (Eq. IT6)) which is used for the first time 
to describe the CP data in various systems. 

Fits to the I- V data at different temperatures as well as 
disorder are shown in Figs. 2,4 and 6, and to the scaled 
curves in Figs. 3,5,7,8 and 9. The near perfect agree- 
ment between the experimental data and theoretical fits 
is quite remarkable in view of range of data covered - 
more than nine orders of magnitude in conductance and 
nearly five orders of magnitude in field. Furthermore, 
fits to the scaled curves using the scaled version of the 
GM-expression (Eq. IT51) provide absolute justification 
of the assumption (Eq. 1171) and underline the relations 
that exist among the coefficients in the GM expression 
but were not foreseen in the theory^. The fact that 
coefficients c„'s used in fits to the scaled curves corre- 
spond closely (within errors) to c„'s (Table II) obtained 
from fittings of individual curves also demonstrate those 
relations in a self-consistent manner. It may be recalled 
that two assumptions about the coefficients, namely Eqs. 
(fT7|) and ([2Tj) led to the important expression Eq. (|22|) 
for the non-linearity exponents xt in the GM-model, 
that allowed a number of conclusions to be drawn on 
the properties of xt itself. Fitted values of the coeffi- 
cients <t„'s (Eq. I16p for different values of n given in 
Table II are plotted against a using log-log scales in 
Fig. 12. The linearity of the plots amply validate the 
assumption (Eq. 1211) . Interestingly, the latter, under 
certain condition, can be derived within the original con- 
text of GM-model^. According to the theory, coefficients 
£„ ~ exp[— 2w/a(n + 1)] are exponentially increasing al- 
beit slowly varying functions of n where w is the typical 
barrier thickness. The Ohmic conductance S is sum 
of two terms - ~ exp(— 2w/a) from direct tunneling 
and Ei ~ exp(— w/a) from resonant tunneling. For thick 
barrier, S D ~ exp(— w/a) so that S n ~ Y, Vn with 



for n > 2. Thus, y n = 0.67, 0.5, 0.4 ... does satisfy 
inequalities for xt > (Eq. 123)1 . In fact, the value of y 2 
i.e. 0.67 is close to some experimental values in Table III 
but seems to be much slowly decreasing than required. 
More importantly, there is a spectrum of values rather 
than a single value for a given channel. 

A close look at Table II reveals yet another intrigu- 
ing feature of channel selection - only the 'even' channels 
(i.e., channels with even number of localized states) ap- 
pear in the fittings to systems like PEDOT, PPy film 
etc. having positive nonlinearity exponents (xt > 0) 
whereas only the 'odd' channels appear in the fittings to 
other systems like PPY (powder) having negative non- 
linearity exponent (xt < 0). As a result, the lowest chan- 
nel contributing to non-Ohmic conduction in the former 
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FIG. 12. (Color online) Log-log plots of GM coefficients vs. 
linear conductances. Solid lines are linear fits to the data with 
slopes as indicated. 



is a two- impurity channel (n a = 2), and a three- impurity 
channel (n Q = 3) in the latter. There was no instance 
of mixing of terms of the two series (i.e. even and odd) 
in any of the samples investigated here. Interestingly, no 
such selection of channels was detected in a model system 
of metal-amorphous silicon-metal tunnel junctions where 
all channels from n—2 to 5 were found to be present in 
some Tj-V curves^. The theory, at present, does not pro- 
vide any clue to this phenomena of selection, in general, 
and to the properties of CP systems responsible for this, 
in particular. Similar selection has been also noticed in 
other hopping systems like carbon nanofibers and amor- 
phous/doped semiconductors^. We also note an oppo- 
site trend for PPy films with different quenched disorder 
(fig. 6b) having a negative nonlinearity exponent where 
the best fits are obtained with even channels in the GM 
expression. This is not surprising as one of the curves 
was earlier part of the data collapse of field dependent 
conductances at various temperatures. The non-linearity 
exponent was found positive (inset, fig. 6a) and fits to 
the GM-expression consisted of even channels only. Dif- 
ferences between the systems described by even and odd 
channels can be seen in the values of c n 's in Table II. 
In case of even channels, C2's are close to 1 and rest of 
c n 's are very small. But in case of odd channels, c„'s are 
comparable as in PPy (powder). 

It is seen from Figs. 5,7 and 8 that the CP's like 
PEDOT, PPy and PDA were described by the GM ex- 
pression with only even non-Ohmic channels: oja = 
1 + c 2 q 4/3 + c 4 g 18 / 5 + c 6 <7 40 / 7 where q = F/F Q . The data 
ranges were sufficiently large that each system needed at 
least two non-Ohmic terms for fitting. While the coeffi- 
cient C2 was close to 1 in all the systems, C4 was 0.011, 
0.0024 and 0.000123 respectively. Such decreasing values 
of c 4 means decreasing rate of rise of conductivity with 
field. One indicator of this trend could be 54, the value of 
q in a system above which the fitting (the dashed curve 
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in the corresponding figure) with only the lowest channel 
(n=2) deviates from the scaling curve. The increasing 
value of (74, roughly 2,10 and 20 respectively in the three 
systems do correlate with C4. More significantly, the de- 
creasing C4 correlates well with the increasing nonlinear- 
ity exponent which are 0.16, 0.23 and 0.51 respectively 
(Table III) : smaller the value of a nonlinearity exponent 
is, the steeper is the conductivity curve (or larger the 
value of C4). This correlation, in fact, is a result of the 
relation Eq. ([7]). The higher order coefficients like c 6 are 
expected to exhibit similar correlation in values. If one 
adopts the view that each distinct value of the nonlin- 
earity exponent constitutes a different universality class 
then one appreciates the unique role that the values of 
c n 's play in giving identity to the corresponding scaling 
function. Incidentally, 54 is given by <r/a = 1 + C294 P2 
from Eq. (1181) . may be interpreted to be the onset 
field for the 4th channel. 

It is quite remarkable that only three nonlinear terms 
were required to describe data in PPy (film) over nine 
decades. The finite number of terms is consistent with 
the requirement that the GM series must terminate at 
some finite n (the point iv in section II). The highest term 
with n=6 (p n — 5.7) was necessary to fit the data at high 
field (F/F a > 1) i.e. data at T ~ 2.1 K which is nearly 
given by a power-law (Fig. 6a). Actually a log-log plot of 
the data yielded a slope of zt = 5.57 close to 5.7. How- 
ever this value of zt is higher than 1/xt = 1/0.23 ss 4.39 
in apparent disagreement with Eq. ([7]) for xt > 0. The 
same happens also in PDA where zt(~ 3.6) > 1/xt = 2. 
Reasons for this quantitative discrepancies are not clear 
although qualitative compliance is obvious as discussed 
above. The relation (Eq. [7]) has been routinely observed 
to hold in amorphous semiconductors^. It is to be noted 
that in both PPy (film) and PDA, F Q at low tempera- 
tures lie above the linear fits (Fig. 10). Whether this is 
indicative of existence of two slopes or not merits further 
careful measurements at low temperatures to resolve this 
issue. 



C. Non-linearity exponents xm 

Non-linearity exponents are the concrete outcome of 
the adopted scaling procedure and are displayed in Ta- 
ble III for various CP systems under different conditions. 
An exponent can be obtained in two ways. One method 
has been already illustrated in Figs. 3, 5, 7, 8 and 9. This 
involves collapsing I-V curves gathered at different val- 
ues of some parameter M (e.g. temperature) by suitable 
choices of F Q and a . The exponent is then obtained us- 
ing Eq. ([5]). This method does not require any knowledge 
of the scaling function and is solely based upon the scal- 
ing property of the I- V curves. The other method makes 
use of the explicit functional form of a I-V curve (or, 
scaling function), when available. It has been shown in 
the above that GM-expressions (Eqs. [TBI and [Tg)) describe 
the relevant data in an excellent manner. Therefore, we 
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FIG. 13. (Color online) Measurement paths along with as- 
sociated nonlinearity exponents in various doped polypyr- 
role samples in Temperature-Disorder plane. Disorder here 
is parameterized by the Ohmic resistivity at 20 K. The data 
marked with star are from Ref. [2r3 . Note that data from a 
sample simultaneously belongs to two intersecting lines. 



can use the expression (Eq. I22[) for obtaining xm in terms 
of model parameters, both direct and derived. y n 's are 
derived from plots like those shown in Fig. 12. There are 
as many y„'s as the number of inelastic tunneling chan- 
nels. According to Eq. (|2"2"1) . each one of them should 
yield the same xt as an yet another test of consistency 
in applicability of the GM-model to CP's. This seems to 
be well borne out within errors (Table III) when the con- 
tribution from the second channel is large enough to yield 
reliable values as in PPy (powder) and PPy (film) (this 
work). Considering the fact that digitization errors are 
not accounted in the values quoted in case of samples of 
other works, agreement between values of the exponents 
obtained using two methods is quite reasonable. In case 
of measurements at different magnetic fields, since con- 
ductivities are related by some constant factors, we have 
y n = 1 for all n. Eq. (|2"2"]l leads to xb = in agreement 
with the value from scaling consideration as discussed 
earlier. It may be noted that y n 's in Table III for both 
positive and negative exponents satisfy inequalities of Eq. 
(|2"3")l . Also, all determined xm's lie between the bounds 
-1/2 and 3/4 (see the point i). Incidentally, putting Eq. 
(gSt into Eq. (J22J) yields x M = l/(n + 2) instead of a 
constant value. 

Clearly, there is no universal exponent for CP's. There 
is no universal exponent even for a given conducting poly- 
mer. As seen in Table III, the same PPy sample exhibits 
two different values of the exponent, 0.23 and depend- 
ing upon the variable M. Three different PPy samples 
exhibits as many values of the exponent xt (-0.33, 0.23 
and -0.17), both positive and negative. This plethora of 
values could simply indicate that xm depends upon the 
path of measurement in the variable space that also in- 
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TABLE III. Comparison of nonlinearity exponents x m obtained using scaling and extended GM analysis (Eq. 1221) in various 
CP systems. M stands for the variable used to vary a : T-Temperature, D-Disorder and B-Magnetic field. y n is given by 



a ° 



with values of a„'s of some samples displayed in Table II. 
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XM 


XM 
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scaling 
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-0.31 ±0.03 




1.78 ± 0.08 
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2.67 ±0.13 
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1.32 ± 0.04 






PPy (film) 
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0.158 ±0.003 


0.170 ±0.005 


0.774 ± 0.007 




~ 0.12 


~ 0.57 


PDA (crystal)*! 


T 


0.51 ±0.02 


0.49 ±0.02 


0.34 ± 0.02 
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0.33 ± 0.02 


0.29 ± 0.06 


0.62 ± 0.08 










FIG. 14. (Color online) Schematic diagram of possible hop- 
ping paths in real space between electrodes for fields (F) less 
(upper) and more (lower), than an onset field (F ). The solid 
circles represent localized states. The lower path corresponds 
to two-step tunneling, the lowest non-Ohmic channel accord- 
ing to Glatzman-Matveev (Ref. Eih . See text for further 
details. 



eludes quenched disorder. Fig. 13 illustrates four such 
paths of measurements in PPy samples along with the 
corresponding exponents in Temperature-Disorder plane. 
Obviously, unraveling of the details of such dependency 
in the variable space has to await theoretical understand- 
ing of xm that is lacking at present. In fact, different 
values of xm in the same system prove that the hopping 
network is in both physical as well as energy space in con- 
trast to the conduction network in composites at room 
temperature^, that is supposed to be purely geometrical 
in real space. The field scale can not be simply deter- 
mined by geometrical topology of disorder alone. Oth- 
erwise, temperature that only provides an energy scale 
could not have an effect on the scale. This is consistent 
with result, xb = since magnetic fields do not change 
either energy levels of charge carriers or disorder, and 
therefore, do not warrant a change in the field scale. 



D. Final picture and Issues with the GM- model 



It is evident that the field-dependent data in CP's are 
in excellent agreement with the GM-expression (Eq. [T6|) 
and the modified one (Eq. [TBI . However, in spite of this 
unexpected success the applicability of the GM-model re- 
mains a moot issue. Let us consider the broad picture of 
transport in lightly doped CP's, that emerges from the 
results presented here and its analysis. The phenomena 
of scaling and the GM-expression as the scaling func- 
tion appears to hold good irrespective of the zero-bias 
conduction mechanism. Except for PA nanofibres, all 
other systems considered in this work exhibit VRH-type 
transport at low bias. In contrast, the model consider- 
ing conduction across a barrier of width w predicts that 
the temperature dependence of conductance in the limit 
eV <C ksT is given by a temperature-independent term 
representing direct tunneling plus a series of terms that 
are sequentially triggered as the temperature is increased. 
The series is similar to Eq. (|16|) with ksT replacing eV 
with the lowest one being proportional to T 4 / 3 . Con- 
sider now the problem of non-Ohmic conduction at a low 
temperature (eV 3> fc^T). In the model, the conduc- 
tion across a barrier proceeds from direct tunneling to 
resonant tunneling to directed hopping along quasi-one 
dimensional paths of localized states and finally, in the 
bulk limit (w 3> Ru) to variable range hopping^.. Thus, 
the VRH that appears in the original GM-model as a 
limiting process appears in the present problem at the 
very outset of low bias. So, how does one envisage the 
evolution of a percolative^ trajectory of a charge carrier 
in VRH (Ohmic) regime as shown schematically in the 
upper panel in Fig. 14 to a linear trajectory in the non- 
Ohmic regime as shown in the lower panel of Fig. 14? 
As the applied field is increased beyond the onset value, 
the conduction becomes non-Ohmic. According to the 
GM-picture, non-Ohmic part in conductance results from 
progressive triggering of the inelastic multi-step tunnel- 
ing channels starting from the lowest two-step tunneling 
(at least when the latter has a non-zero contribution) 
shown in the lower panel of the figure. The conductance 
corresponding to two-step tunneling is proportional to 
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exp(— w/3a) where w is the barrier width. Obviously, it 
can not be a macroscopic length (w 3> a) typical of the 
samples considered here. If we take w to be of the or- 
der of disordered region along a polymer fibril, then one 
has to consider many such processes throughout a sam- 
ple in comparison of only one in the original model. If 
VRH applies to the whole sample it is not clear what will 
be the total trajectory corresponding to two-step tunnel- 
ing. Furthermore, since Ohmic-non-Ohmic transition is 
a continuous process it is hard to imagine a transforma- 
tion of paths of a charge carrier as indicated in the figure. 
To summarize, one has to reckon with the fact that the 
expression (Eq. [16]) that was supposed to describe hop- 
ping across a thin film only is now required to describe 
hopping across apparently multitudes of such 'films' em- 
bedded inside bulk systems. 

Apart from the conceptual problems outlined above, 
we list here the interesting features involving particu- 
larly the coefficients a n that were observed and need to 
be explained theoretically. Firstly, the incidence of chan- 
nel selection. This is not a universal phenomenon. It 
occurs in CP's and probably, other materials but is ab- 
sent in tunnel junctions^. It is necessary to understand 
the exact conditions in the GM-model that give allow 
such a selection phenomenon, and then how those condi- 
tions may be satisfied in CP's. Secondly, the assumption 
(Eq. [TT|) that makes scaling analysis possible within the 
GM-model. Thirdly, Eq. (12T1) that relates each channel 
coefficient to the linear conductivity. Simple arguments 
yield only y n < 1 for n > 1 (Eq. [2HJ). But y n > 1 is 
necessary to have negative nonlinearity exponents. 

VI. CONCLUSION 

In this paper, we have reported field-dependent con- 
ductivities in various CP systems as a function of tem- 



perature and quenched disorder. We demonstrated that 
each of various conducting polymer systems possess a 
single field scale and exhibits the associated scaling. A 
phenomenological scaling equation that led to extraction 
of nonlinear exponents was used to analyze the nonlinear 
transport data. It was argued that experimental evidence 
points to the linear conductivity as a legitimate scaling 
variable. Surprisingly, the GM-expression for multi-step 
tunneling proves to be an excellent fit to the I- V curves 
as well as the scaled curves. A couple of assumptions 
are made to make the GM-model compatible with scal- 
ing. Experimental values of the exponents fall within 
the predicted limits of -0.5 and 0.75. A theory capable of 
explaining the nonlinear exponents, particularly negative 
ones, is lacking. The value of the exponent depends upon 
the path of measurement in the variable space. Some is- 
sues concerning applicability of the GM-model to CP's 
have been discussed. 

The scaling that has been observed in CP's has been 
also found in many other disordered systems including 
composites. It is believed that such scaling may be truly 
an universal feature of disordered systems particularly 
with localized states. All the samples considered here 
are three-dimensional. It will be interesting to know how 
such scaling fare in lower dimensions. 
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